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ABSTRACT 


New  second-  and  third-order  algorithms  are  presented  for  cal¬ 
culating  translating  and  rotating  steady-state  solutions  of  the  2D 
Incompressible  Euler  equations  (which  we  call  V-states).  These  are 
piecewise  constant  regions  of  vortlclty  and  the  contours  bounding 
them  are  obtained  by  solving  Iteratively  a  nonlinear  Integro-dlfferentlal 
equation.  New  limiting  contours  with  corners  are  obtained  and  compared 
with  local  analytical  solutions.  The  precise  results  correct  mistakes 
for  limiting  contours  that  were  previously  given. 


STEADY-STATE  SOLUTIONS  OF  THE  EULER  EQUATIONS  IN  TWO  DIMENSIONS; 

ROTATING  AND  TRANSLATING  V-STATES  WITH  LIMITING  CASES 
I.  NUMERICAL  ALGORITHMS  AND  RESULTS 

» 

H.,  M.  Wu  ,  E.  A.  Overman,  II,  and  N.  J.  Zabusky 

Institute  for  Computational  Mathematics  and  Applications 
Department  of  Mathematics  and  Statistics 
University  of  Pittsburgh 
Pittsburgh,  PA  15261 

1.  INTRODUCTION 

Solutions  of  the  two-dimensional  Incompressible  Euler  equations 
will  elucidate  properties  of  very  large  Reynold’s  number  flows,  as 
may  occur  In  planetary  atmospheres  and  oceans  The  method  of  contour 
dynamics.  Introduced  by  Zabusky,  Hughes  and  Roberts  [1],  provides  a 
computationally  convenient  approach  because  a  two-dimensional  problem 
Is  reduced  to  one-dlmenslon.  That  Is,  In  contour  dynamics  the  sources 
of  the  flow  are  plecewlse-constant  regions  of  vortlclty,  ^Ich  we  call 
PAVR’s  (finite  area  vortex  regions)  or,  equivalently,  the  contours 
bounding  these  regions. 

The  parameters  that  describe  the  range  of  existence  and  stability 
of  steady-state  FAVR  configurations  may  elucidate  properties  of  these 
flows.  Klrchoff  [2]  found  that  an  elliptical  FAVR  was  a  steady-state 
solution  and  Love  [31  Investigated  Its  stability.  No  other  closed 
form  solutions  have  been  found.  Deem  and  Zabusky  [4]  found  new  steady- 
state  FAVR  configurations  by  solving  numerically  a  nonlinear  integro- 

*Permanent  Address:  C«q>uter  Center,  Academia  Slnlca,  Beijing,  China. 
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differential  equation  for  the  contours  using  a  Newton-Raphson  pro¬ 
cedure.  For  example,  they  found  several  Isolated  rotating  states  of 
m-fold  symmetry  that  are  bifurcations  from  harmonic  waves  on  a  cir¬ 
cular  FAVR  and  one  Isolated  symmetrically-shaped  dipolar  (l.e., 
oppositely-signed  vortlclty)  translating  state.  Examples  are  sketched 
In  Fig.  1.  They  considered  these  regions  of  plecewlse-constant 
vortlclty  as  members  of  larger  seta  of  steady-state  solutions  and  referred 
to  them  generally  as  *’V-states".  Saffman  and  colleagues  also  applied 
this  technique  to  calculate  shapes  of  doubly-connected  rotating  V-states 
[5]  and  periodic  states  modeling  free  shear  layers  [6]  and  wakes  [71. 
Plerrehumbert  [8]  applied  an  efficient  first-order  relaxation  method 
and  obtained  twelve  members  of  the  set  of  symmetrically-shaped  dipolar 
translating  states.  Plerrehumbert  and  Wldnall  [91  also  applied  this 
algorithm  to  calculate  free  shear  layer  models.  Burbea  and  Landau  [10] 
applied  the  same  algorithm  and  obtained  further  examples  of  m-fold 
symmetric  rotating  V-states  for  3  £  m  ^  6.  In  both  [8]  and  [10]  the 
limiting  V-states,  where  the  contour  Is  nonanalytlc,  are  defective. 

This  occurs  because  their  algorithms  seem  unable  to  handle  singular  and 
near-singular  contours  and  because  the  spatial  resolution  Is  Inadequate 
In  regions  of  large  curvature.  In  the  latter  paper,  this  results  In 
errors  In  the  nvimerlcal  calculation  of  the  parameter  range  of  existence 
of  rotating  V-states.  The  analytical  calculation  by  Burbea  [11]  of 
this  range  Is  also  Incorrect  and  both  calculatloris  are  discussed  In 
Section  6. 

I 

In  this  paper  we  present  a  new,  fast,  accurate  and  computationally 
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the  neighborhood  of  nonanalytlcity .  For  the  translating  V>state, 
there  are  two  possible  limiting  cases:  the  two  regions  may  touch 
at  one  point  or  they  may  have  a  common  boundary.  Analytically,  we 
have  been  able  to  exclude  the  former  case.  For  the  latter  case  we 
have  established  analytically  that  the  two  (one-sided)  tangent  angles 
at  a  nonanalytlcal  point  may  differ  only  by  ir/2  (l.e.,  a  corner). 
Numerically,  we  have  confirmed  the  existence  of  this  solution  as 
shown  In  Fig.  4b.  For  the  rotating  V-states  we  have  established 
analytically  that  the  tangent  angles  at  a  nonanalytlcal  point  may 
differ  only  by  0  (l.e.,  the  tangent  angle  Is  contlnuouat)  or  ir/2. 

Numerically,  we  have  confirmed  the  existence  of  the  ir/2  corner  for 
3  5  m  ^  6  as  shown  In  Pig.  7. 

In  a  recent  letter,  Saffman  and  Tanveer  [12]  also  did  a  local 
analysis  of  the  limiting  translating  case  and  obtained  a  ir/2  corner. 
They  also  used  numerical  methods  to  calculate  this  state,  but  only 
provide  a  gross  figure  and  Insufficient  Information  to  allow  a  detailed 
comparison  of  results. 

In  Section  2  we  present  analytical  preliminaries  and  derive  two 
nonlinear  Integro-dlfferentlal  equations  for  the  boundary  which  are 
the  basis  for  our  new  second-  and  third-order  accurate  algorithms.  In 
Section  3  we  analyze  the  limiting  cases  In  regions  where  a  contour  may 
become  non-analytlcal  and  prove  the  claims  made  above.  In  Section  4 
we  present  the  discretized  versions  of  the  Integro-dlfferentlal  equations 
and  iterative  algorithms  for  obtaining  both  translating  and  rotating  V- 
Btates.  In  Sections  5  and  6  we  discuss  properties  of  the  numerical 
solutions  for  translating  and  rotating  V-states,  respectively.  In  both 
cases,  magnified  views  of  the  contours  are  given  In  the  region  of 
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nonanalytlcal  behavior  and  they  are  compared  to  the  local  solutions 
of  Section  3.  We  also  present  a  thorough  discussion  of  the  sensitivity 
of  the  results  In  the  neighborhood  of  nonanalytlcal  points  to; 
algorithms;  discretization  procedures;  error  criteria;  and  boundary 
conditions . 

2.  THE  INTEORO-DIPPERENTIAL  EQUATIONS  FOR  V-STATES 

The  Euler  equations  can  be  written  In  vortlclty-stream  function 
form  as 

+  uw^  +  vw„  »  0,  (2,1a) 

w  X  y 

where 

(2.1b) 

and 

(u,v)  2  (2.1c) 

If  the  vortlclty  Is  represented  by  a  set  of  N  piecewise  constant 

c 

functions  of  strength  In  regions  Dj  with  boundaries  we 

can  express  the  stream  function  as 

♦  (x»y)  ■  I  w*  ff  0(x-C,y-.n)dCdTi,  (2.2)/ 

J-1 

irtiere  Q  Is  the  Green's  funetlMi  for  the  Laplaelan  In  the  unbounded 
dcMUtln 


0(x-€,y-n)  --(2»)“Hogt(x-€)^  ♦  --<2ir)“hog  L, 

(2.3) 
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If  Green's  theorem  is  applied  to  the  result  of  substituting  (2.2) 

Into  (2.1c)  we  obtain  an  expression  for  the  velocity  as  a  sum  over 
the  N  contour  Integrals,  namely 

w 

Nc 

(u,v)  =  (u(x,y)  ,v(x,y))  «  (2Tr)“^  I  [u],  [  log  ^(dC,dn),  (2.^1) 

J-1  ^ 

where  [wjj  la  the  Jump  in  vorticlty  (outslde-lnslde)  at  3Dj  and 
where  the  dependence  on  time  has  been  suppressed.  If  we  Integrate 
by  parts  we  obtain 

Nc 

(u,v)  -  (2ir)“^  I  [m],  f  ^“^(x-e,y-n)d^.  (2.5) 

J-1  ^ 

The  contours  ai'e  assumed  to  be  piecewise  Llapounov,  where  a  Llapounov 
curve  Is  one  which  possesses  a  unique  continuous  tangent  angle,  a, 
but  not  necessarily  a  curvature,  at  each  point  [13, 14].  Thus,  we 
require  that  each  contour  consists  of  a  finite  number  of  segpnents, 
each  of  which  possess  a  unique  a  at  every  point  but  may  have  an 
Infinite  curvature  at  the  ends. 

Kelvin's  theorem  requires  that  a  particle  on  the  boundeu*y  remains 
on  the  bovuidary.  Hence,  for  steady>state  solutions 

~  *  ^particle  "  -  *  -boundary*  (2.6) 

where  Kpjixtlcle  "  ®  outward  normal  to  the  con¬ 

tour.  For  translating  V-states  (2.6)  can  be  written  as 

u  sin  a  -  (.v-V)oo8  a  «  0,  (x,y)  e  30^  and  dDg, 


(2.7T) 
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where  V  is  the  translational  speed  (the  velocity  is  in  the  y-direction 
from  Pig.  la).  For  rotating  V-states  (2.6)  can  be  written  as 

u  sin  a  -  V  cos  a  »  -fiR  dR/ds,  (x,y)  e  (2.7R) 

where  n  Is  the  angular  velocity  of  the  state,  s  la  the  arclength, 
and  R(s)  Is  the  radius  from  the  origin  to  the  contour  as  shown 
In  Pig.  lb.  Note  that  for  simplicity  in  (2.7)  and  henceforth  we 
label  equations  with  '•T'*  (translating)  or  "R"  (rotating)  according 
to  the  state  being  considered.  We  now  assume  that  the  contours  are 
"star shaped",  where  the  single-valued  R(0)  is  defined  with  respect 
to  a  convenient  origin  as  shown  in  Pig.  1,  so  that 

(x,y)  •  (R(0)cos  0  +  R(0)sln  0), 

where  (2.8) 

R(0)  •  [(x(0)  -  Xq)^  +  y^(0)]^^^ 

(Xq  ■  0  for  the  rotating  case).  Hence,  (2.7)  becomes 

u  dy/d0  -  (v-V)dx/d0  -  0,  (x,y)  e  or  aDg,  (2.9T) 

or 

u  dy/d0  -  V  dx/de  +  12R  dR/d0  ■  0,  (x,y)  e  3D.  (2.9R) 

Also,  we  can  use  (2.8)  to  write  (2.9)  explicitly  in  terms  of  R(0) 

as 

dR/de  -  XR, 


(2.10) 
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where 


or 


X  »  [(V-v)sln  0  -  u  cos  6]/[(V-v)cos  9  +  u  sin  0]  (2.11T) 


X  »  -[u  cos  0  +  V  sin  0]/[u  sin  0  -  v  cos  0  +  OR].  (2. HR) 


As  discussed  In  Section  4  we  will  use  (2.9)  for  our  second-order 
scheme  and  (2.10)  for  our  third-order  scheme. 

An  alternative  "stream  function"  form  of  the  Integral  equations 
Is  obtained  by  writing  (2.9)  as 


Oy'l'Xdy/de)  +  (V  +  3^ij<)(dx/d0)  »  0,  (x,y)  e  or  aDg,  (2.12T) 

or 

(3  i|/)(dy/d0)  +  O  i(/)(dx/d0)  +  flR  dR/dO  »  0,  r  »  R(0),  (r,0)  e  3D, 

y  * 

(2.12R) 

and  Integrating  to  obtain 

i(»(x,y)  +  Vx  ■  Cj,  (x,y)  e  3Dj,  J  «  1  or  2,  (2.13T) 

or 

i|»(r,0)  +  (Q/2)R^(0)  -  c,  r  -  R(0),  (r,9)  e  3D.  (2.13R) 


Plerrehumbert  [8]  used  (2.13T)  and  Burbea  and  Landau  CkS  used  (2.13R) 
to  obtain  steady-state  solutions,  whereas  we  use  the  "velocity"  form 
as  described  above. 

In  this  paper  we  solve  two  classes  of  problems: 
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PROBLEM  T:  (symmetric  dipolar  translating  V-state) 

This  state  Is  symmetric  about  both  axes  as  shown  In  Fig.  la  with 
vortlcity  0)^  ■  +1  and  *  -!•  Then,  given  Xg  =  1  and  0  £  x^  £  1 
find  R(9)  and  V  that  satisfy  (2.9T)  or  (2.10)  with  (2.11T). 

PROBLEM  R;  (m-fold  symmetric  rotating  V-state) 

This  state  with  vortlcity  +1  has  m  Identical  sef  rs  each  of 
which  has  two  reflectlonally  symmetric  subsectors  as  she  '  j.n  Fig.  lb 
for  m  »  3.  Then,  given  m,  Rg  -  R(ir/m  -  ir)  =1  and  R^  =  R(-ir/2)  >  1 
find  R(0)  and  0  that  satisfy  (2.9R)  or  (2.10)  with  (2. HR). 
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3.  ANALYSIS  FOR  LIMITING  V-STATES 

3.1  Introduction 

In  this  section  we  summarize  the  local  analysis  of  V-states,  pre¬ 
sented  In  detail  elsewhere  [15].  We  use  a  local  expansion  to  obtain 
the  equation  for  the  boundary  of  a  symmetric  dipolar  translating  or 
m-fold  symmetric  rotating  V-state  In  the  neighborhood  of  a  possible 
singularity.  The  analysis  provides  a  necessary  condition  for  the 
behavior  of  the  boundary.  Namely,  the  difference  In  tangent  angles 
at  a  singular  point  can  be  only  ir/2  (l.e.,  a  corner)  or,  for  the 
rotating  case,  0  (l.e.,  the  tangent  angle  Is  continuous).  For  the 

corner  V-states  we  obtain  the  local  equations  of  the  boundary  and 
compare  them  with  the  numerical  results  in  Sections  5  and  6.  In 
Section  3.4  the  analysis  Is  applied  to  Klrchoff's  elliptical  vortex 
to  validate  the  procedure. 

The  general  method  of  analysis  Is  as  follows .  For  convenience  we 
use  polar  coordinates,  (r,0),  with  the  singularity  at  the  origin .  As 
shown  In  Fig.  2  we  assume  that  the  V-state  can  be  oriented  so  that 
It  Is  symmetric  about  the  y-axls  and  lies  in  the  upper  half  plane 
(l.e.,  the  vortlclty  on  the  negative  y-axls  Is  0).  We  do  all  cal¬ 
culations  In  the  right  half-plane,  l.e.,  -ir/2  £  6  £  x/2,  using  symmetry 

to  complete  the  V-state.  We  assume  In  some  neighborhood  of  the  origin, 

0  <  r  <  6,  that  the  boundary,  0  ■  0(r),  is  once  continuously  dlfferentlabl 
so  that  the  tangent  angle  Is  continuous.  (However,  this  does  not  pre¬ 
clude  11m  d9(r)/dr  ■  +  «.)  We  expand  the  integral  expression  for 
r-*-0  “  2 

i|»  In  (2.2)  In  terms  of  r  to  0{r  log  r).  Since  the  V-state  Is 
stationary  In  the  appropriate  reference  frame,  the  value  of  the  stream 
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function,  t()®.  Is  constant  on  Its  boundary  and  so,  without  loss  of 
generality,  V»®(r,8(r))  -  0  for  0  <  r  <  6.  We  write  9  (r)  as 

0(r)  -  Oq  +  ej^(r)  (3.1) 

where  Bq  is  the  tangent  angle  of  the  (right  half  of  the)  V-state 

at  the  origin,  l.e.,  11m  0(r)  ■  0^,  and  so 

IM-O  “ 

11m  0, (r)  »  0.  (3.2) 

r^O  ^ 

Since  the  solution  of  Poisson’s  equation,  (2.1b),  Is  once  con¬ 

tinuously  differentiable  in  all  of  3R^,'  we  can  expand  it  In  a  Taylor 
series  with  remainder,  or 

•s 

0  -  i<)®(r,9(r))  -  i|/®(r,0Q+03^(r))  «  t|»®(r,9Q)  +  3Qi»»®(r,a)0j^(r) ,  (3.3) 

where  a  =  o(r)  Is  In  the  open  interval  between  0q  and  0q  +  0^(r). 
Thus, 

03^(r)  -  V(r,0o)/3gt;^®(r,a(r)).  (3.^») 

We  obtain  the  possible  values  of  0q  from  (3.2),  l.e., 

11m  i|»®(r,0Q)/3g»|>*(r,a(r))  -  0,  (3.5) 

r^O  ” 

where  0g  can  be  substituted  for  a(r)  In  (3.5)  since  11m  a(r)  ■  0«. 
In  those  cases  where  ♦  (r,0g)  $  0,  which  Includes  the  comer  cases. 
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> 


s  2 

we  then  expand  ip  up  to  0(r  )  and  solve  for  0j^(r)  from  (3.^). 

(To  lowest  order  In  r  we  can  replace  a(r)  by  Oq.) 


3.2  Translating  V-states 

The  limiting  symmetric  translating  V-state  consists  of  two  PAVR’s 
with  a  common  boundary  on  the  y-axls  and  with  vortlclty  ■  +1  and 
0)2  ■  -1.  The  right  FAVR  Is  composed  of  two  contour  segments,  namely 
9(r)  as  described  above  and  0  «  ir/2  which  Is  the  common  boundary. 

At  the  end  of  this  subsection  we  will  show  that  the  two  PAVR’s  of  a 
llmltlng-case  V-state  cannot  touch  at  only  one  point.  (However,  we 
cannot  rule  out  the  possibility  of  an  Isolated  V-state  whose  PAVR's 
have  one  point  In  common.) 

For  0  <  6  «  1  and  r  <<  6,  we  find  that  [15] 


(r,0)  -  [V 


"  11 


COS  $  dpd<f]r  cos  0 


-  (Mir)  ^(1  +  cos  20f.)r‘^log  r  sin  20  +  o(r‘^log  r). 


(3.6) 


where  V  Is  the  speed  of  the  V-state.  Note  that  '('®(r,0)  ■  0(r) 

If  the  leading  term  of  (3.6)  Is  non-zero.  Thus,  from  (3.5)  cot  0q  «  0 
so  Qq  "  i  leading  term  vanishes,  l.e.. 


V 


cos  ^  dpd^. 


(3.7) 


then  (3.5)  yields  tan  20q  -  0  so  0q  ■  0  or  +  ir/2.  We  first 
discuss  0Q  *  0. 

.For  Sq  •  0,  from  [15] 
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(|;®(r,0)  »  -  (2ir)“^r^log  r  sin  26  +  (2ir)“^(C^+C2)r^  sin  20 


-  C^ir)“^r^  20  cos  20  -  sin  26  + 


+  o(r  0^(r)) , 


,  if  it/2  >  6  >  0 


IT  COS  20,  if  0  >  0  ^  -Tt/2 


(3.8) 


where 


II 

D,(6) 


sin  2«  dpd^  +  (1/2) (1  +  21og  6), 


(3.9a) 


Q 

^2  "  “f  0j^(p)dp, 


(3.9b) 


and  D^(6)  Is  that  subset  of  whose  distance  from  the  origin  Is 

>  6.  Prom  (3.^)  we  find  that 

0^(r)  -  (■ir/4)/(-log  r  +  +  C2)  +  o((log  r)"^).  (3.10) 

Substituting  (3.10)  Into  (3.9b)  we  obtain  a  quadratic  equation  for 
C2  and  find  that 


C2  -  %  [log  6  -  +  (dog  6  -  -  (irV4»^]+  0((log 


The  slngvilar  behavior  of  the  curve  near  the  origin,  where 


(3.11) 


6^(r)  «  1  can  be  seen  from  the  slope. 
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^  «  sin  9(r)  r  cos  6(r)  6'(r)  ^  .  2 

tx  cos  6\.r)  -  r  sin  e(r)  6*(r)  “  +  O(0^(r)), 


(3.12a) 


and  the  curvature. 


K  M  2e«(r)  +  r^(9»(r))3  +  r9»(r) 
(1  +  r^(0'(r))2)3/2 


(8/ir)9j(r)/r  +  O(0^(r)/r). 


(3.12b) 


Note  that  the  tangent  angle  is  0(l/iog  r)  and  the  curvature  Is 
"1  —2 

0(r  (log  r)  )  so  that  llm  ic  »  <».  Also,  the  velocity  at  any 

x^O 

point  (r,9)  is 


-  r  log  r  cos  9  +  ir”^  (Cj^+C2)r  cos  0 


-  (2v)  r  -20  sin  9 


{V  sin  0,  v/2  >  0  >  0  I" 
-V  sin  0,  0  >  0  >  -v/2/ 


V  -  -  ir~^  r  log  r  sin  0  +  (C^+C2)r  sin  0 


+  o(r  0, (r)). 


(3.13a) 


-  (2ir)"^r[20  cos  0  +  IT  cos  0]  +  o(r^  0j^(r)), 


(3.13b) 


so  that  on  the  curve  u  -  0(r  log  r)  and  v  ■  0(r) 
For  Bq  "  Cl5]  that 


♦*(r,0)  -  Cq  r  CO8  0  +  (2ir)*^  aln  20  +  o(r^  cos  0), 


(3.14) 
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where 


and 


-Ml 


Cq  ■  V  +  IT  *  ]  ]  cos  ♦  dpd^ 


■  II  p“^  aln  2^  dpd^. 


(3.15a) 


(3.15b) 


Substituting  6  ■  0Q  +  0^(r)  In  (3.14)  we  obtain 

0  -  i|(®(r,eQ+ej^(r) )  «  -Cq  r  sin  e^(r)  -  (2»)“^ 

+  o(r^  sin  6^(r))  . 


sin  20j^(r) 

(3.16) 


If  we  divide  (3.l6)  by  r  sin  02(r)  and  let  r  0,  we  find  that 
Cq  ■  0.  If  we  then  divide  (3.16)  by  sin  B^(r)  and  let  r  0 
we  find  that  ■  0.  However,  »  0  In  (3.15b)  If  and  only  of 
the  origin  of  the  coordinate  system  Is  on  the  horizontal  line  through 
the  centroid  of  .  That . Is ,  the  two  PAVR ' s  touch  at  only  one  point . 
This  follows  because  the  Integrand  of  (3.15b),  Is  euitl symmetric 

about  ^  ■  0  and  the  V-state  Is  assumed  to  be  symmetric  about  this 
horizontal  line  of  symmetry  of  D^.  If  we  do  the  same  analysis  for 
0Q  ■  -ir/2  we  obtain  the  same  result. 

We  now  show  that  this  configuration  leads  to  a  contradiction. 

If  the  analysis  Is  repeated  with  the  origin  on  this  line  of  symmetry 
[13],  then  we  again  obtain  (3.14)  and  so,  again,  Cq  •  0  which  contra¬ 
dicts  the  value  obtained  from  (3.15a).  To  see  this  we  use  the  numerical 
solutions  obtained  in  Section  5  (or  In  [8]).  We  estimate  the  Integral 
in  (3.15s)  by  taking  the  llBd.t  as  0  of  a  circle  of  radius 
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(l-x^)/2  centered  at  ( (l+x^)/2,0) .  Prom  Pig.  4  this  circle, 
which  we  denote  by  C(x^)',  lies  Inside  the  corresponding  V~state, 
denoted  by  Thus, 


II 


cos  ^  dpd^ 


=  11 


COS  ^  dpd(fr  ■  (Tr/2)(1  - 


x^)V(l  +  x^). 


(3.17) 


and  so 


II 


cos  ^  dpd^  >  1/2. 


(3.18) 


Prom  Table  I,  0  >  V  >  -.258  for  all  the  V-states  and  so  Cq  >  0.2. 
(Actually  Cq  »  0.6  for  x^^  ■  10"^.)  Thus,  solutions  with  PAVR’s 
touching  at  one  point  are  excluded  so  the  only  possible  solution  Is 
the  comer  solution  shown  In  Pig.  3a. 


3.3  Rotating  V-states 

For  the  rotating  V-state  we  let  ■  Ug  ■  +l  In  Pig.  2  and  find 
In  [15] 

\t»®(r,e)  ■  II  sin  ♦  dpd^  -  ny]r  sin  0 

Dl 

+  (4w)”^  sin  20q  r^  log  r  cos  20  +  o(r^  log  r),  (3.19) 

i^ere  0  Is  the  angular  speed  of  the  V-state  about  the  centroid  at 
(0,y).  Note  that  i|*®(r,e)  ■  0(r)  if  the  leading  term  in  (3.19)  Is 
non-zero.  Thus,  from  (3.5)  tan  6q  ■  0  so  0q  ■  0.  If  the  leading 
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term  vanishes,  l.e.. 


II 


sin  ^  dpd^. 


(3.20) 


then  (3.5)  yields  cot  29q  *  0  so  0q  “  +  ir/^.  We  will  Investigate 
only  Sq  -  +  ir/4,  since  the  numerical  results  In  Section  6  Indicate 
that  0Q  ■  -ir/4  Is  not  a  limiting  case.  The  two  solutions  0q  “  0 
and  ir/4  are  shown  In  Pig.  3b. 

For  0^  ■  ir/4,  from  [15] 


i|»°(r,0)  ■  (4Tr)~^r^  log  r  cos  20  +  (2ir)“^CCj^+C2)r^cos  20 


ry[- 


2fffl  +  (l/2)cos  20  +  (ir/2)0  gin  20 


,  If  ir/2 
Sin  20,  If  tt/4 


^  0  >  ir/4  1  " 
>  0  >  -v/z/ 


+  o(r  0^(r)), 


(3.21) 


where 


’l  "  \\  "  (1/^)(1  +  21og  6), 

I),(6) 


(3.22a) 


5 


sin*  6^(p)dp. 


(3.22b) 


Thus,  from  (S.'t) 


0^(r)  -  CO-(3/8)]ir/(log  r  +  2(C^+C2))  +  o((iog  r)“*). 


(3.23) 
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Substituting  (3.23)  Into  (3.22)  we  obtain  a  quadratic  equation  for 
Cg  and  find  that 

C2  -  -(1/4)  [log  6  +  2C^  +  ((log  6  +  2C^)^  -  8it^ [n-( 3/8 ) ]^ )**] 

+  0((log  «)“^).  (3.24) 

In  this  case  the  slope  Is 

dy/dx  -  1  +  2ej^(r)  +  0(e^(r)),  (3.25a) 

while  the  curvature  Is 

K  .  -2([R-(3/8)]ir)‘^e2(r)/r  +  0(e3ir)/r),  (3.25b) 

which  Is  similar  to  the  translating  case,  Eqs.  (3.12).  The  velocity 
Is  also  similar  to  that  previously  given. 

3.4  Elliptical  V-state 

The  analysis  In  the  rotating  case  above  Is  also  valid  for 
analytical  V-states  since  the  only  constraints  we  have  used  are  that 
the  origin  lies  on  the  V-state  and  the  V-state  Is  symmetric  about  the 
y-axls.  We  will  now  consider  the  rotating  case  when  >  0  to  show 
that  we  obtain  the  correct  result  for  ellipses. 
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r(r,e) 


r  sin  6  +  (2*  )"^(Cj^^)r^  cos  26 


-  (l/4)r^ 


-2fl  + 


{: 


-(1/2)co8  20,  ir/2  >  0  >  0 


(l/2)cos  28 


,  ir/2  >  0  >  0^1  - 

>  +  0(r3), 

,  •  0  >  0  >  -ir/^ 


where 


Cq  -  ir“^  II  sin  (|)  dpd(f-ny. 


(3.26) 


(3.27a) 


-  Jl.  . 


cos  2(^  dpdi^. 


D^(6) 


and 


6 

'  Cg  -  -(1/2)  I  p’^  sin  203^(p)dp. 
0 

Prom  (3.^)  we  obtain 


(3.27b) 


(3.27c) 


0j^(r)  ■-(2Cq)  ^  (Cj^  +  Cg  +  mi  -  (ir/4))r  +  O(r^).  (3.28) 

Substituting  this  Into  (3.27c)  we  obtain 

°2  "  (2Cq)~^(c.,  +  irn  -  ir/4)6  +  0(6^),  (3.29) 

so 

6j(r)  -- (2^0^^  (Cj^  +  irfl  -  Tr/4)(1  +  (2CQ)“^fi)r  +  O(rfi^).  (3.30) 

From  Section  6  we  find  that  for  an  analytical  V-state  Cq  >  0  and 
Cq  ♦  0  as  the  limiting  V-state  Is  approached  so  that  <-►  +  «•.  At 
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the  limiting  V-state  the  equation  of  the  curve  Junps  to  (3.23). 

We  now  show  that  If  the  equation  for  aui  ellipse  In  our  coordinate 
system. 


_  _  2sln  0  / /cos^  6  L  sln^  e^  ^ 

r  *  - -  t  2 (3.31) 

/a  b 

Is  expanded  near  the  origin,  the  result,  r  s  2a^e/b,  agrees  with 
(3.30).  First  we  substitute  (3.31)  Into  (3.27a)  and  (3.27b)  and 
obtain 

Cq  “  ifb  ~  (3.32a) 

and 

°1  "  j(t)/a^)8f(|  ir)[|^  -  1-3  +  (3.32b) 

For  the  ellipse  »  ab/(a+b)^  [2]  and  If  we  substitute  for  n,  C,. 
and  In  (3.30),  we  obtain  0j^(r)  ■  (|- b/a^)r  +  0(r5^)  which 

agrees  with  the  expansion  of  (3.31)  for  0  <<  1.  (Note  that  If 
a  ■  b  then,  for  arbitrary  n,  0j^(r)  ■  r/a  +  0(r6^)  as  It 

should  since  the  circular  FAVR  Is  a  V-state  for  any  n.) 

3.5  Summary 

In  this  section  we  have  examined  the  behavior  of  a  V-state  In 
the  neighborhood  of  a  singularity.  To  apply  the  analysis  we  put 
the  origin  of  the  coordinate  system  at  the  point  In  question  and  then 
require  that  the  V-state  can  be  oriented  to  be  symmetric  about  the 
y-axls  as  In  Fig.  2.  This  can  be  done  for  both  the  sytmnetrlc  dipolar 
translating  V-state  and  tha  m~fold  symmetric  rotating  V-state  as  shown 
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In  Figs.  4b  and  7.  (The  analysis  can  also  be  done  without  assuming 
S3nnmetry.)  In  the  translating  case  we  find  that  the  difference  In 
tangent  angles  at  a  singularity  can  be  only  ir/2.  We  examine  this 
comer  case  further  and  find  that  for  r  <<  1  the  curvature,  k. 

Is  0(r“  (log  r)”  )  which  Indicates  the  difficulty  In  numerically 
calculating  the  V-state.  In  the  rotating  case  we  find  that  the 
difference  In  tangent  angles  at  a  singularity  can  be  0  or  ir/2. 
Again,  for  the  corner  case  k  *>  0(r  (log  r)"  ).  Prom  numerical 
results  0  corresponds  to  an  analytical  V-state  and  ■n/2  to  the 
limiting  state  for  3  <  m  <  6  (and  we  assume  for  all  m) . 
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4.1  Velocities 

We  present 
They  both  use  a 
on  the  contours 

-  (2> 

where 

and  the  mean  pos 

and 

4.2  Second-Orde 

The  second- 
latlng  V-states 
V-etates  for  1. 
comer.  It  oonv 
algorltha  to  be 
Halt  Inc  nMianal 


4 .  NUMERICAL  ALGORITHMS 


ew  second-  and  third-order  accurate  algorithms, 
econd-order  accurate  representation  of  the  velocities 
(u,v)jjj,  namely  a  trapezoidal  discretization  of  (2.5) 


Nc 

-1  j 

J-1 


I 

1«1 


(J) 


k_M+V 


(y^-h 


UT735 


1+15 


(/(J) 

'‘^k,l+l  ^k,l^* 


(4.1) 


)  e  is  number  of  points  on  3D 

tlons  are  given  by 


*  Cyk-ni^')^]’*. 


Algorithm 

rder  accurate  algorithm  Is  used  to  obtain  trans- 
—7 

or  10  '  jc  <  0.90  and  rotating  m-fold  symmetric 

§  § 

5  £  <  R^(m),  where  is  the  value  at  the 

rgea  mov9  rapidly  than  the  third-order  Iterative 
resented  in  the  next  subsection.  However,  for  the 
tleal  contours  the  third-order  algorithm  ^Ives  more 
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accurate  results  as  shown  In  Sections  5  and  6. 

We  use  the  .second-order  accurate  discretization  of  (2.9), 

“ic+ii‘yk  -  ■  °>  1  i  k  1 H. 

or 

“k+tj^^k  -  ''ktls^'^'k  *  '“^2’  1  1  1  N. 

Where  and  defined  in  (4.2). 

Here,  N  +  1  is  the  number  of  points  on  the  segment  of  the  contour 
for  -ir/2  £  6  £  0  in  the  translating  case  and  -it/2  £  6  £  v/ra  -  ^ 
in  the  rotating  case  (the  solid  lines  in  Pig.  1).  We  define 

*  R(ej^)(cos  ej^,sin  0j^)  (4.5) 

where  ■  0  for  the  limiting  translating  V-state  and  all  rotating 
V-states  and  *q  *  j  smalytical  translating  V-states. 

Substituting  (4.5)  into  (4.4),  we  obtain 

Rj^’ -  Pk+jjRit+i  •  0.  1  <  k  <  N,  (4.6a) 

or,  alternatively, 

\  -  PkiijRic.i  -  0.  2  <  k  <  N+1,  (4.6b) 


(4.4T) 


(4.4r) 


where  Pjj+jj  Is  defined  as 


,  .  “k*».  °kti  - 

k+>i  -  sin  9^  - 


'ikijLl’fflVl 


(4.7T) 


or 
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p 


- 


Sin  cos  6^^^^  ^ 

sin  cos  +  (n/2)Rj^ 


C^.7R) 


To  obtain  convergent  algorithms,  we  find  it  necessary  to  use  a 
three-point  scheme  and  a  relaxation  procedure.  First  we  average  (4.6a) 
and  (4.6b)  to  obtain 

-  ? ’'kUVl  *  \  -  I  SilciN.  (4.8 


where 


and  “  Rg*  This  discrete  representation  of  our 


nonlinear  integro-differential  equation  can  be  solved  for  R^^  if  we 
know  Thus,  if  we  have  Just  completed,  the  n-th  iteration  we 

know  R^”^  and  so  can  find  ^“k+V^k+%^  ('^.D.  We  then  cal¬ 

culate  the  new  velocity  by  summing  (4.4)  to  obtain 


or 

Thus,  we  can  calculate  by  (^-7)  and  solve  the  linear  equation 

(4.8)  for  Rj^  -*>  We  obtain  by  "relaxing"  R^"^  and 

by 

-  W*  *  (1-P»)R^"^  (4.10) 

where  u*  ■  0.6.  We  discuss  the  initial  guess  used  In  the  appropriate 
section. 


4.3  Third-Order  Algorithm 

Chronologically,  we  first  obtained  the  second-order  algorithm 


given  in  Section  4.2  but  found  that  In  the  limiting  cases  It  could 
give  Inadequate  results.  A  third-order  algorithm  Is  readily  obtained 
by  using  the  differential  equation  for  R,  Eq.  (2.10),  and  weighting 
three  adjacent  Rj^  In  the  manner  described  below.  First,  we  dis¬ 
cretize  (2.11)  using  a  midpoint  method 

ARj^/A0k  - 

We  rearrange  and  obtain 

\  "  °k+A+l  “  1  <  k  <  N, 

or 

Here  is  defined  as 

°k+i5  -  -  I  ^k+is^®k^'^^^  I  ^k+55^®k^»  (4.13a) 

where 

^k+3a  =  ®^k+35*  ®^k+35^ 

(4.13b) 

and  X  is  defined  in  (2.11).  The  third-order  property  is  achieved 
by  weighting  R^  In  the  following  manner: 

*  "He  -  -  O’  2 1  “  i  <“•»  > 

The  local  analysis,  carried  out  in  Appendix  A,  shows  that  terms 
O((A0k)^)  cancel  exactly  if 

-$k(A0u.i)^  +  (l-Bk)(A0k)^  -  0. 


(4.12a) 


(4.12b) 


(4.15) 
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(The  algorithm  discussed  In  Section  4.2  could  be  made  third-order 
by  the  same  type  of  procedure  but  It  was  easier  to  expand 
(4.13),  rather  than  (4.7).) 

In  the  translating  case  It  was  found  that  convergence  could  not 
be  obtained  for  the  limiting  V-state,  even  for  the  third-order 
algorithm,  when  the  boundary  condition  at  the  singularity  was  a  fixed 
angle  as  opposed  to  a  fixed  point ♦  (This  Is  discussed  further  In 
Sections  5  and  6.)  The  solution  oscillated  over  a  small  range  In 
the  neighborhood  of  the  singularity. 

Hence,  to  obtain  convergence  we  use  a  two-step  procedure:  a 
method  of  stabilization  [l6]  followed  by  a  method  of  relaxation.  For 
the  method  of  stabilization  we  replace  LRj^  *  0  In  (4. l4)  by  a  dis¬ 
cretized  version  of  “  0  where  ^  sis  t  •*■  <*>.  That 

Is,  we  solve  for  the  (n+l)st  Iteration  by 

L(R^"^)  +  y(R^^'^^^  -  R^"^)  =  0,  (4.16) 

where  we  have  Introduced  the  Intermediate  variable 
y  ■  0.1.  We  again  obtain  "relaxing"  R^^^ 

by  (4.10).  (Note  that  (4.10)  and  (4.l6)  are  readily  combined  Into 
one  equation  In  out  program.)  If  ■  1/2  In  (4.l4)  the  algorithm 
Is  second-order  accurate  but  not  Identical  to  the  algorithm  In 
Section  4.2.  We  have  not  uded  this  second-order  algorithm  In  this 


paper . 
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4.4  Summary  and  Convergence  Criteria 

The  second-  and  third-order  algorithms  are  summarized  as 


follows : 

(1) 

Compute 

/  (n)  (n)x 
^'^k+35*^k+i5^ 

from  using 

(4.1) 

(2) 

Compute 

or 

from  (4.9). 

(3) 

Compute 

(4.7)  or  from 

(4.13)  for  the  second 

or  third-order  algorithms,  respectively. 


(4)  Compute  from  (4.8)  for  the  second-order  algorithm  or 

^(n+l)  (4.16)  and  from  (4.10)  for  the  third-order 

algorithm.  In  either  case  use  Gaussian  elimination  to  Invert 
the  trldlagonal  matrix. 

(5)  Continue  the  iteration  until  the  error  criterion  Is  satisfied. 

A  run  Is  terminated  If 


N+1 

- 

k-1 


(n+1) 

k 


<  £ 


(4.17) 


where 


e 


^5  10“’^ 
5  X  10”^ 


for  translating  states 
for  rotating  states. 


V. 

VThen  convergence  Is  obtained  we  find  that  the  original  Integro 

-8 

differential  equations  are  satisfied  to  5  x  10~  ,  l.e.. 


I  ‘  5  >■ 

or 

I'^k+Js^^k  "  ’^k+55^*k  *  (fi/2)AR^|  <  5  ^  10“®. 

This  accuracy  was  verified  on  the  DEC-10  (a  36-blt  machine) 
at  the  University  of  Pittsburgh  by  continuing  runs  In  double 
precision  once  the  required  accuracy  was  obtained  In  single- 
precision. 
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5.  NUMERICAL  CALCULATION  OP  TRANSLATING  V-STATES 

5.1  Analytical  V-atates 

To  compute  the  sequence  of  states  with  10  5.  £  0.90,  we 

use  the  second-order  algorithm.  For  case  1  the  Initial  state  Is  a 
half-circle  of  radius  0.05  centered  at  (Xq,0)  ■  (0.95,0),  with  N  +  1 
equally  spaced  nodes,  l.e.,  AO  ■  ir/N.  For  the  remaining  cases  In  the 
sequence,  the  Initial  state  Is  obtained  by  expanding  linearly  the 
previously  obtained  solution  with  (Xq,0)  »  (^aCl+x^)  ,0) .  For  cases 
1-14,  when  the  error,  e,  (4.17),  £  10~^  we  adjust  the  nodes  so 

that  the  distance  between  the  adjacent  nodes  Is  Inversely  proportional 
to  the  curvature,  *  '^k+Js*  makes  the  local  error  the  same 

In  each  Interval  [17].  For  cases  15-17,  A0j^  Is  the  same  as  obtained 
for  case  l4.  With  this  discretization,  we  continue  Iterating  until 
convergence  Is  obtained. 

The  results  obtained  with  the  second-order  algorithm  and  N  ■  120 
are  summarized  In  Figs.  4  and  5  and  Table  I.  Fig.  4a  represents  one 
sector,  l.e.,  1/4  of  the  V-state,  for  cases  1  (x^  -  0.90)  through 

_ii 

13  (x^  ■  10  ).  To  show  the  power  of  the  algorithm  we  have  enlarged 

the  scale  by  »  400  In  Pig.  4b  and  show  cases  l4  (x^  -  10”^)  through 

^*7 

17  (x^  “  10  ').  We  observe  that  the  contours  are  nested  and  tend  to  a 
limiting  contour,  the  lowest  In  Fig.  4b,  discussed  below. 

In  Fig.  5  and  Table  I  we  present  properties  of  the  sequence  of 
states  where  A  Is  the  area  of  one  side,  P  Is  the  perimeter,  x 
Is  the  X  coordinate  of  the  center  of  area,  V  Is  the  translational 
speed,  V/Vq  Is  the  normalized  translational  spded  where  Vq  ■  A(4it5E)"^ 
Is  the  translational  velocity  of  two  point  vortices  with  circulation 
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+  A  and  separation  2x,  R  •  (A/ir)^,  H  ■  maxjyl  Is  the  maximum 
vertical  extent  of  a  sector,  and  a  ■  2H/(l-x^)  Is  the  aspect  ratio. 

The  dots  In  Pig.  5  are  the  results  of  Plerrehumbert  [8]  and  the  com¬ 
parison  Is  excellent  except  for  H  In  the  limiting  case^as  discussed 
below.  The  convergence  criterion  that  e  <  Is  obtained  with 

the  second-order  algorithm  In  less  than  70  Iterations.  An  Iteration 
step  with  N  »  120  requires  8  seconds  of  CPU  on  the  DEC-10.  Most  of 
this  time  Is  consumed  calculating  the  velocities  at  the  nodes.  A  thorough 
discussion  of  accuracy  and  sensitivity  Is  given  In  the  following  subsectlor 

5.2  Limiting  V-state  (x^  «  0) 

As  Indicated  In  Pig.  4b  the  V-states  tend  to  a  limiting  state.  In 
Section  3  we  observed  that  a  limiting  contour  could  approach  the  y-axls 
only  when  the  tangent  angle  at  the  axis,  a^.  Is  0.  In  this  subsection  we 
Investigate  the  sensitivity  of  this  approach  angle  with  the  second-  and 
third-order  algorithms.  The  following  paragraphs  discuss  the  boundary 
conditions.  Initialization  and  discretization  of  this  nonanalytlcal  state. 

We  assume  that  the  boundaries  of  both  contours  of  the  limiting  V- 
state  lie  on  the  y-axls  from  (0,-y»)  to  (0,y*).  We  let  the  center 
of  our  polar  coordinate  system  be  at  (0,0)  and  compute  the  velocities 
(ujj,Vjc)  In  two  parts.  Plrst,  we  do  a  niunerlcal  Integration,  Eq.  (4.1), 
for  -ir/2  1  0  1  0  as  previously  and,  second,  we  do  an  analytical  In¬ 
tegration  of  Eq.  (2.4)  from  (0,0)  to  (0,-y*)  as  discussed  In  [l]. 

Two  types  of  boundary  conditions  are  used  at  9  ■  -ir/2.  Plrst, 
to  find  the  comer  solution  we  set  j  e*_ir/2  "  ®  fitting  a 

quadratic  polynosilal,  symmetric  about  the  y-axls,  through  the  second 
and  third  points.  The  resulting  matrix  can  be  transformed 
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to  a  tridlagonal  form  and  solved  as  previously.  This  limiting  case 

Is  shown  In  Pig.  ^b.  We  call  this  the  comer  boundary  condition. 

Second,  we  fix  the  point  on  the  y-axls,  R(-ir/2),  which  enables  us 

to  determine  the  sensitivity  of  this  limiting  V-state. 

Because  of  Its  singular  character,  as  discussed  In  Section  3,  the 

limiting  case  Is  approached  very  slowly  and  the  selection  of  a  ”good" 

Initial  state  is  Important.  As  described  in  Appendix  B  the  initial 

state  is  derived  by  smoothing  the  last  analytical  state  (i.e..  No.  17, 

—7 

x^  *  10  ' )  with  a  high  density  of  nodes  near  the  corner  as  shown  in 
Pig.  4b.  We  find,  using  the  third  order  algorithm  and  the  corner  boundary 
condition,  that  y*  ■  1.66855,  1.66898  and  I.6691I  for  N  ■  30,  60  and 
120,  respectively.  The  convergence  Is  very  slow  as  the  solution 
exhibits  a  damped  oscillation  around  y»  and  requires  2i3»000  Iterations 
to  satisfy  the  error  criterion.  The  second-order  algorithm  solutions  do 
not  converge  but  simply  oscillate  slowly  about  y*  with  a  range  of 
siO.OOl. 

The  calculations  of  y*  are  consistent  with  the  fact  that  the 
algorithm  Is  actually  only  second-order  accurate  near  the  y-axls  because 
(u,v)  cannot  be  expanded  In  a  Taylor  series  at  the  singularity  (see 
(3.13)).  Using  second-order  RlchardAon  extrapolation  on  N  -  60  and 
120,  we  find  to  five  significant  figures 

1.6691  1  y»  <  1.6692.  (5.1) 

Note  that  Plerrehumbert 's  limiting  V-state  has  a  cusp  for  the 
singularity  and  y*  ■  1.705  (■  3.41/2).  It  seems  to  us  that  his  dis¬ 
tribution  of  nodes  In  the  nelghbomood  of  the  singularity  was  In- 
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adequate  since  we  could  obtain  "V-states"  with  ■  -7r/2  emd 
+Tr/2  when  the  neighborhood  of  the  singularity  was  Inadequately  re¬ 
solved.  (In  Pig.  ^b,  Ax^  »  1.4  X  10“^.) 

In  Pig.  6a  we  have  plotted  the  limiting  V-state  In  the  neighbor¬ 
hood  of  the  singularity  for  N  »  60  and  120  (the  dots).  This  Is  a 
magnification  of  «.  30  over  Pig.  4b.  Note  that  In  Pig.  4b  the  tangent 
angle  at  x  »  0  does  not  seem  to  be  0  (even  with  a  magnification  of 
fls.200)  but  Is  seen  to  be  much  closer  to  0  In  Pig.  6a  (a  magnification 
of  «  6,000),  which  shows  the  singular  nature  of  the  curve. 

Also  In  Pig.  6a  we  plot  the  equation  for  the  curve  (the  solid 
line)  from  Section  3,  Eq.  (3.10),  with  parameters  and  C2  obtained 

from  the  solution  with  N  ■  120.  has  been  calculated  numerically 

from  (3.9a)  and  Cg  from  (3.11).  For  6  -  0.000014,  +  Cg  -  0.1196 

while  even  for  6  ■  0.049,  ^2  “  *1106.  In  Pig.  6b  we  continue 

this  comparison  on  a  larger  scale  to  show  the  quality  of  the  asymptotic 
formula.  Also,  the  velocity  of  the  V-state  as  calculated  from  (3.7) 

Is  -0. 25797  while  numerically  It  Is  -0.25793. 

To  determine  the  sensitivity  of  the  eLLgorlthm  we  use  the  second 
boundary  condition,  i.e.,  fix  R(-ir/2),  near  y*.  In  Table  II  we  show 
the  results  using  both  the  second-  and  third-order  algorithms  for 
N  ■  60  where  the  x  coordinates  of  the  nodes  correspond  to  those  In 
Pig.  6a.  Note  that  In  all  the  cases  the  solution  tends  very  rapidly 
to  the  comer  solution.  (We  have  given  7  sl^if leant  figures  for 
comparison  purposes,  but  trust  only  the  first  5.)  For  example,  the 
maximum  difference  In  is  1.673000  -  1.668973  “  0.004027  and  In 
is  0.000060  (where  x^  «  0.00013).  In  all  runs  using  the  fixed- 
point  boundary  condition, o<»iver(enoe  is  obtained  In  slOO  Iterations. 


32 


6.  NUMERICAL  CALCULATION  OP  ROTATING  V-STATES 


Aneilytlcal  V-states 


We  compute  the  sequence  of  states  for  1.05  1  with  the 

second-order  algorithm,  where  R*  Is  the  value  obtained  for  the  limiting 

V-state.  Our  Initial  state  Is,  for  -ir/2  <  0  <  6  -  j  it, 

—  —me 


®  1 

IS  c»|i  9J 


R^°^(0)  -  +  (1-Ra^*  (2-9  ) 


(6.1) 


where  0^  ■  ir/m,  where  9  ■  (6+|-Jf)/0jn  and  A0j^  Is  constant.  Pig.  7 
shows  (x,y)  and  curvature  plots  for  one  analytical  V-state  and  the 
limiting  V-state  for  3  <  m  <  6.  In  the  curvature  plots  the  abscissa 
Is  the  arclength  scaled  so  that  It  Is  0  at  R^  and  1  at  Rg.  The 
properties  of  the  analytical  V-states  are  given  In  Table  Ilia.  The 
limiting  cases  (•)  are  given  for  comparison  and  discussed  below. 


6.2  Limiting  V-state 


e  1 

CAP  6J 


For  the  limiting  V-state  our  Initial  state  Is 


R^°^(0)  -  Rj^  +  (1-Ra^  9^(3-29) 


(6.2) 


where  R^  ■  1.73#  1.32  and  1.24  for  m  ■  3#4,5  and  6,  respectively. 

The  angular  difference  is  either  constant  as  In  the  analytical 

case  or  Increases  nearly  linearly  with  k  as  discussed  In  Appendix  B. 

For  the  latter  we  start  with  A ^  ■  0.1*#  0.01*  or  0.001*.  We  again 
use  two  types  of  boui^ary  conditions  at  6  ■  -ir/2:  the  corner 
boundary  condition,  dr/dxL  ,  boundary 
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condition,  l.e.,  R(-ir/2)  ■  R^.  This  corner  boundary  condition  Is 

obtained  by  using  a  linear  combination  of  either  the  first  3  or  4 

points,  as  described  In  Table  Illb. 

Our  results  for  3  £  m  <  6  are  contained  In  Table  Illb  for  the 

corner  boundary  condition.  Since  the  total  number  of  points  on  the 

V-state  Is  2mN  the  time  required  for  each  iteration  Increases  as 
2 

m  and  so  we  will  only  consider  the  case  m  «  3  In  detail.  At  the 
end  we  will  comment  about  the  other  cases. 

We  perform  a  sensitivity  study  for  m  »  3  using  various 
algorithms,  boundary  conditions  and  discretizations  as  shown  in  Table 
Illb.  With  the  "linearly”  Increasing  discretization  *  0.1*,  0.01“ 

and  0.001“)  we  find 


n*  -  0.30120  +  0.00004, 
R|  -  1.7352  +  0.0003. 


(6.3) 


In  the  rotating  case,  unlike  the  translating  case,  both  the  second- 
and  third-order  algorithms  converge  to  the  limiting  V-state  for  the 
comer  boundary  condition.  The  second-order  algorithm  converges  In 
s  500  Iterations  while  the  third-order  algorithm  requires  ■=.  2000 
Iterations  because  the  stabilization  and  relaxation  procedures  delay 
the  convergence.  Since  the  second-order  algorithm  Is  also  much  faster 
with  the  fixed-point  boundary  condition,  we  use  only  the  third-order 
algorithm  when  high  accviracy  Is  required. 

The  range  of  existence  of  for  2  <  m  £  6  Is  shown  as  the 

solid  vertical  lines  In  Fig.  8  and  fl*,  the  lower  end  of  these  lines. 


can  be  fit  with 
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n.  -  i  (m-2)[^ 


m  -  1.195 


m*"  -  2.071m  +  .2085 


(6.4) 


where  |n*  -  fi*|  <  10“^.  In  previous  analytical  work  [11]  the 
lower  end  of  the  range  of  existence  was  given  as  (1/2) (m-2)/(m-l) , 
shown  In  Pig.  8  as  the  dots.  This  result  Is  incorrect  because  Burbea 
linearized  about  the  circular  V-state  and  Interpreted  his  results  as 
being  valid  in  the  nonlinear  region.  In  previous  computational  work 
DcQ numerical  solutions  were  obtained  below  our  range  of  existence. 

For  example,  for  m  *  3  they  presented  a  solution  at  0  «  0.2822 
which  had  regions  of  negative  curvature.  These  Incorrect  results  are 
probably  due  to  Inadequate  discretization  procedures  and  to  the  fact 
that  spurious  "solutions"  can  be  obtained  for  >  R»,  as  we  will 
discuss  below. 

Since  the  range  of  existence  of  was  missed  previously  [ID] , 

we  present  the  results  of  several  sensitivity  studies  for  m  «  3  and 
R^  near  R^.  First,  the  dots  In  Pig.  9  show  the  maximum  curvature 
as  a  function  of  R^  obtained  with  the  third-order  algorithm,  N  *>  .60 
and  a  discretization  of  0.01®  for  Rg  <  1.70  and  0.001®  for 
Rg  >  1.70.  The  x's  are  due  to  Burbea  and  Landau  Do]  and  we  have  not 
plotted  their  last  value  of  k  ■  236  at  R^  ■  1.923,  which  is  well 
to  the  right  of  the  figure. 

Second,  we  used  the  fixed-point  boundary  condition,  both  algorithms 
and  various  discretizations  to  obtain  the  results  In  Fig.  10a,  where 
we  have  plotted  the  tangent  angle  at  the  singularity,  i.e, ,  vs  R^. 

The  solid  line  shows  the  small  range  obtained  for  R^,  1.735  <  R^  <  1.736, 

using  the  third-order  algorithm  and  0.01®  (0.001®  Is  not  shown  since 
It  Is  undistlngulshable  from  0.01*).  Also  note  that  all  solutions  cross 
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■  45*  at  »  1.735.  In  Pig.  10b  we  show  the  behavior  of  the 

contours  near  R|  using  the  third-order  algorithm  and  0.001“  where 
\  -  1.73333,  1.73493  and  1.73559. 

As  shown  In  Pig.  10a,  we  can  also  obtain  "solutions”  for  R^  >  R^ 
and,  indeed,  for  x  90* .  This  Is  due  to  the  fact  that  the  various 
discretized  forms  of  Euler's  equations  have  different  solutions  than  the 
continuous  equations,  Eqs.  (4.1).  For  larger  values  of  R^  we  still 
can  obtain  convergence,  but  we  find  that  the  solutions  behave  in  a  strange, 
algorithm-dependent,  noncontinuous  fashion.  For  example,  for  the  third- 
order  algorithm  with  N  «  60  and  0.001*  the  solution  Jumps  back  from 
s  90*  to  «.  0*  as  R^  Is  Increased  slightly  (to  si. 7358).  For 

the  remaining  cases,  R(02)  decreases  as  R^  increases  until  for 
R^  ■  1.80  we  find  that  ®  1  2  k  <  N+1  so  that  the 

solution  looks  like  a  circle  with  a  sharp  spike.  We  take  this  as  evidence 
that  we  have  passed  the  range  of  existence  of  steady-state  solutions  to 
the  continuum  equations,  which  do  not  have  V-states  with  cusps. 

Finally,  In  Fig.  11  we  compare  the  numerical  results  (the  dots)  for 
the  0.01*  run  with  the  4  point  boundary  condition  in  Table  Illb  to  the 
formulas  In  Sec.  3,  Eq.  (3.23).  has  been  calculated  numerically 

from  (3.22a)  and  then  Cg  from  (3.24).  For  6  »  0.00044  then 

+  C2  •  -O.692  while  even  for  6  «  0.085,  +  C2  ■  -0.694.  The 

comparison  Is  excellent  for  the  entire  sector,  -■n/2  £  0  £  -ir/6.  Also, 
a  Is  0.30122  from  (3.20)  while  It  Is  0.30121  In  Table  Illb.  In  Fig.  12 
we  compare  an  analytical  V-state,  Eq.  (3.30),  with  the  solution  of 
the  third-order  algorithm,  N  ■  60  and  0.001  for  R^  ■  1.7349.  In 
this  case  the  equation  requires  fl  as  well  as  the  nodes  on  the  con¬ 
tour  from  the  numerical  calculation.  Using  the  value  of  0  <■  0.301170 
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we  find  that  Cq  *  0.00052,  »  -M.4489  and  the  curvature  at 

obtained  by  differentiating  (3-30)  Is  ic  »  1495*  which  agrees  with 

\ 

the  numerical  solution  to  three  significant  digits.  We  consider  all 
of  the  above  a  sufficient  validation  of  our  calculations  of  the  limiting 
V-states . 

The  only  difficulty  we  have  encountered  for  3  £  m  £  6  Is  that 
the  curvature  oscillates  near  0  »  -ir/2  as  the  limiting  V-state  Is 
approached.  In  fact,  the  calculations  of  the  V-states  In  Pig.  7  were 
done  with  the  nodes  equally  spaced  In  angle  (see  Table  Illb)  and  for 
m  »  4,5  and  6  It  Is  Just  possible  to  see  wiggles  In  the  curvature  plots 
near  the  singularity  (l.e.,  s  »  0).  (For  m  »  3  It  Is  possible  to 
remove  the  oscillations  by  a  Judicious  choice  of  discretization  while 
for  4  £  ra  £  6  the  size  and  location  of  the  oscillations  change  with 
the  discretization  but  do  not  disappear.)  However,  for  3  £  m  jc  6,  the 
curvature  does  have  the  correct  sign  at  0  »  -Tr/2.  That  Is,  from  (3.23) 
the  curvature  at  should  be  +  «»  for  m  ■  3  and  4  since  (3/8)  -  n  >  0, 

while  It  should  be  -  *  for  m  >  5  since  (3/8)  -  JJ  <  0.  With 
sufficient  nodes  we  believe  these  oscillations  would  disappear. 
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7.  DISCUSSION  AND  CONCLUSIONS 

We  have  presented  new  accurate  (and  fast)  algorithms  and  refined 
procedures  for  computing  ssnnmetrlc  translating  aind  rotating  V-states 
of  the  Euler  equations  In  two  dimensions.  These  Include  limiting  non- 
analytlcal  contours  with  corner  singularities  that  are  compared  with 
analytical  solutions  [15].  The  agreement  Is  excellent!  These 
singularities  were  missed  In  previous  numerical  work  [8,10], 

Burbea  and  Landau  (10]  proposed  calling  the  rotating  V-states 
"nonlinear  Kelvin  waves".  However,  "V-states"  (vortex  states)  seems 
more  appropriate  since  there  are  already  at  least  two  types  of  Kelvin 
waves  aud  since  Deem  and  Zabusky  first  showed  their  existence  [4]  and 
coined  the  expression. 

In  all  but  the  limiting  cases  the  second-order  accurate  algorithm 

converges  to  the  V-state  quickly,  both  In  the  number  of  iterations 

2 

felOO)  and  the  CPU  time.  The  CPU  time  per  Iteration  Is  0(N  )  compared 

to  Newt on-Raph son  which  is  0(N  ).  The  algorithms  of  Plerrehumbert  [8] 

2 

and  Burbea  and  Landau  [10]  are  also  0(N  ). 

Our  development  of  the  third-order  algorithm  may  seem  ad  hoc  but 
came  about  In  a  search  for  an  accurate  method  to  calculate  limiting 
V-states.  We  have  used  refined  procedures  to  validate  our  results  In¬ 
cluding  various  discretizations  in  the  neighborhood  of  the  singularity 
and  two  boundary  conditions.  In  this  paper  we  have  not  attempted  to 
find  procedures  to  minimize  computation  time.  However,  since  most  of 
the  time  In  an  Iteration  Is  taken  up  calculating  the  velocities.  In 
recent  work  we  have  recalculated  them  every  20  Iterations  when 
£  <  10  .  We  find  It  reduces  the  computation  time  by  a  factor  of  3. • 


•  • 
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We  are  presently  using  the  new  algorithms  and  have  obtained 
asymmetric  translating  V-states  and  V-states  with  nested  contours. 
The  latter  Involves  the  solution  of  coupled  integro-differentlal 
equations,  one  for  each  contour.  However,  there  Is  a  constraint 
that  the  velocities  V  or  fi,  for  all  contours  must  be  equal. 
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APPENDIX  A;  DERIVATION  OP  A  THIRD  ORDER  ALGORITHM 


Since  the  velocities  calculated  by  the  trapezoidal 

rule  (4.1),  they  are  accurate  to  second  order.  Here,  we  show  that 
we  can  obtain  a  third-order  algorithm  by  properly  weighting  two 
equations. 

Let  X  and  R  be  the  solutions  of  R*(0)  *  XR,  (2.10),  (where 
primes  denote  differentiation  with  respect  to  6),  and  Rj^  = 
be  the  solution  of  our  self-consistent  discrete  representation  (4.11). 

Let 

^k+Js  “  ®k+35  ’ 

where  *  O((A0j^)  ).  Then  after  some  algebra  (4.12)  yields 

<A-2) 

where  Gjt+jj  -  expand  Rj^  in  (4.11)  about  Rj^  we 

obtain 

K  -  °k«.*k+i  -  *  0({Ae^)S,  (A.3). 

where 

E(0)  -  (R”»/24)  -  (XR"/8).  (A. 4) 

If  we  substitute  (A. 2)  into  (4.11a)  and  subtract  (A.3)  we  Obtain 

CStA’  -  W"k«-Sk«,>  ♦  w«k  “k«,Vi,  *  ®k«,'''«k>^  ■  0<<iek>''>- 

(A. 5) 


We  apply  the  same  technique  to  (4.11b)  and  find 


<v\>  -  -  ° 


(A. 6) 


To  remove  the  leading  order  error »  O((A0)^),  In  (A. 5)  and  (A.6),  we 

multiply  (A. 4)  by  1  -  and  (A. 5)  by  to  obtain 

■®k°k-is^’’k-l"\-l*  ■  *^”®k'°k+!!*’*k+l"\+l' 

*  [Sk<-ek-!,‘®k-l  <A-1  *  ^ 

«ktA«  *  '^k+>s<''«k>^)3  =  Oide)").  (A.7) 

4 

Thus,  we  must  choose  6j^  so  that  the  term  In  brackets  Is  0((Ae)  ), 
namely 

-6j^(A0j^^l)^  +  (l-Bj^)(A0j^)^  -  0.  (A. 8) 

^  «T_  ^ 

This  follows  because,  to  lowest  order,  •  1>  Rjj_l  ■  \+i* 

=k->s  ■  °<''®k-l>^  «k+l5  •  °<''®k>^  (for  some  c). 

where  the  latter  three  expressions  are  valid  If  the  contour  Is  analytic 


In  the  region  from  k  -  1  to  k  +  1. 


ST 


»  « 


APPENDIX  B:  INITIAL  DATA  AND  DISTRIBUTION  OF  THE  NODES  FOR  LIMITING 
CASES . 

B.l  Translating  Case 

In  order  to  obtain  the  Initial  approximation  R^^^S)  we  begin 
with  the  nodes,  I  ^  1  ^  1  N+1},  from  the  17th  state,  l.e., 

_7 

f  »  10  We  find  the  value  of  k  at  which  yj^  Is  a  maximum,  say 

k  «  K.  Then  we  modify  all  ^  1  by 

*k  “  Xjj^3C<1c-1)/(K+2)]^*5,  1  <  k  <  K+2, 

to  obtain  a  nearly  geometric  ratio  and 

^k  *  1  1  k  1  K+2, 

where  P(x)  Is  the  unique  quadratic  function  satisfying  7^+3  “ 

^K+S  "  *  °- 

B.2  Rotating  Case 

The  Interval  -ir/2  <  0  <  S^-ir/2,  where  S  ■  ir/m.  Is  divided  Into 

—  —  m  m 

N  -f  1  angles  by 

®k  “  N+1. 

For  m  -  3  and  N  -  60  If  A0j^  -  0.01“  then  L  -  691,  A0jj  -  6.2“ 
and  the  ratio  of  the  largest  AO  to  the  smallest  Is  s  620.  If 

AB^  ■  0.001“,  then  L  ■  9950,  A6jj  ■  8.5*  and  the  ratio  Is  a  85OO. 


Properties  of  the  Translating  V-states 


»  *  «  • 
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TABLE  II 


The  Limiting  Translating  V-state  In  the  Neighborhood  of  the 

Corner  (N  ■  60). 


Algorithm 

“1 

«1 

«2 

«3 

R4 

89.3“ 

1.665000 

1.668862 

1.668942 

1.668964 

37.6“ 

1.669001 

1.669015 

1.669012 

1.669010 

2nd 

-0.8“ 

1.669015 

1.669014 

1.669012 

1.669009 

-13.2“ 

1.669021 

1.669014 

1.669012 

1.669009 

-89.7" 

1.673000 

1.669162 

1.669069 

1 

1.669040 

89.7* 

1 

1.664974 

1 

1.668814 

1.668894 

1.668915 

42.0“ 

1.668956  ' 

1.668973 

1.668970 

1.668967 

3rd 

-0.5“ 

1.668973 

1.668972 

1.668970 

1.668967 

-23.4“ 

1.668980 

1.668972 

1.668970 

1.668967 

i 

-89.7“ 

1.673000 

1.669162 

1.669069 

1.669040 

1.668972 

1.669006 

1.669005 

1.669005 

1.669023 

1.66892^ 

1.668963 

1.668^963 

1.668963 

1.669023 


•  • 


Table  Ilia 


Propertlea  ot  Rotatlng^  V-statea 


m 

"a 

Si 

A 

P 

1/Ra 

3 

1.05000 

0.33310 

3.2967 

6.4440 

0.95238 

1.22128 

0.32843 

3.7967 

7.0463 

0.81881 

1.39256 

0.31996 

4.2382 

7.7082 

0.71810 

1.56384 

0.30948 

4.5909 

8.4021 

0,63945 

a 

1.73512 

0.30122 

4.7291 

9.0652 

0.57633 

4 

1.05000 

0.37460 

3.2947 

6.4488 

0.95238 

1.14709 

0.37153 

3.5684 

6.8154 

0.87177 

1.24418 

0.36620 

3.8046 

7.2310 

0.80374 

1.34127 

0.35938 

3.9870 

7.6811 

0.74556 

» 

1.43836 

0.35395 

4.0546 

8.1380 

0.69524 

5 

1.05000 

0.39945 

3.2927 

6.4555 

0.95238 

1.11587 

0.39705 

3.4726 

6.7259 

0.89616 

1.18174 

0.39313 

3.6263 

7.0389 

0.84621  ' 

1.24761 

0.38815 

3.7426 

7.3859 

0.80153 

» 

1.31348 

0.38425 

3.7842 

7.7534 

0.76134 

6 

1.05000 

0.41596 

3.2909 

6.4641 

0.95238 

1.09855 

0.41395 

3.4193 

6.6827 

0.91029  I 

1.14711 

0.41087 

3.5284 

6.9388 

0.87176 

1.19566 

0.40701 

3.6097 

7.2271 

0.83636 

• 

1.24421 

0.40407 

3.6375 

7.5415 

0.80372 

' 
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TABLE  mb 


Properties  of  the  Limiting  Rotating  V-states 


N 

algorithm 

BC 

discretization^ 

SI* 

m 

R» 

m 

A 

1/R* 

m 

60 

2nd 

3 

1*  (C) 

0.30126 

1.7331  ' 

4.7273 

0.57701 

60 

2nd 

3 

0.1"  (L) 

0.30123 

1.73^9 

4.7280 

0.57640 

60 

1 

2nd 

3 

0.01“  (L) 

0.30122 

1.7349 

4.7281 

0.57639 

60 

' 

3rd 

4 

0.01“  (L) 

0.30121 

1.7353 

4.7308 

0.57626 

60 

3rd 

4 

0.001“(L) 

0.30117 

1.7355 

4.7337 

0.57620 

120 

2nd 

3 

t  (L) 

0.30124 

1.7349 

4.7283 

0.57640 

60 

2nd 

3 

0.75*  (C) 

0.35392 

1.4381 

4.0541 

0.69537 

60 

2hd 

3 

o 

O 

« 

O 

0.35397 

1.438:3 

4.0551 

0.69528 

60 

2nd 

3 

o 

• 

ON 

o 

O 

0.38419 

1.3137 

3.7840 

0.76123 

60 

2nd 

3 

0.01“  (L) 

0.38429 

1.3133 

3.7842 

0.76143 

60 

2nd 

3 

0.5*  (C) 

0.40399 

1.2446 

3.6376 

0.80350 

60 

2nd 

3 

0.01“  (L) 

* 

0.40411 

1.2441 

3.6375 

0.80382 

C  ■  Conetant,  L  ■  "Linearly"  increasing 

Corresponds  to  0.01*  with  N  -  60  and  with  an  additional  node  midway 
between  these  nodes. 


FIGURE  CAPTIONS 


Schematic  and  notation  for  V-states  (the  lines  of  symmetry 
are  dotted,  (a)  Translating.  Dipolar  vortlclty,  <^2 
(b)  Rotating.  The  FAVR  shovm  has  3-fold  symmetry. 

Schematic  and  notation  for  a  representative  FAVR  used  In 
the  analysis  In  Section  3.  (Only  In  this  section  Is  the 
origin  of  the  coordinate  system  at  the  singularity.) 

Schematic  showing  the  local  behavior  of  the  limiting  V-state 
near  the  singularity  at  the  origin.  The  possible  values  of 
0Q  are  described  In  Section  3.  (a)  Translating.  The 

vertical  line  Is  the  line  of  Intersection  of  the  two  PAVR's. 
(b)  Rotating. 

(a)  A  sector  of  the  translating  V-states  for  cases  1  through 
13  given  In  Table  I.  (b)  A  magnified  view  of  the  V-states 
In  the  region  of  high  curvature  for  cases  14  through  17  and 
the  limiting  case  (the  lowest  curve).  The  dots  are  the  nodes 
used  In  the  numerical  calculation. 

Global  properties  for  one  contour  of  the  translating  V-state: 
A  (area);  P  (perimeter);  x  (x  coordinate  of  the  centroid); 
and,  V/Vq  (the  normalized  speed).  The. dots  are  the  corres¬ 
ponding  values  from  [8]. 


A  magnified  view  of  the  corner  for  the  limiting  trans¬ 
lating  V-state.  (a)  The  dots  are  the  numerical  solutions 
obtained  with  the  3rd  order  algorithm  and  the  corner  boundary 
condition  for  N  ■  6o  and  120.  The  solid  line  Is  the  local 
solution,  (3.10),  fit  to  the  N  »  120  solution  (b)  The 
comparison  of  Figure  6a  on  a  larger  scale. 

Rotating  V-states  and  their  curvatures,  ic(s),  for  3  <  m  £  6. 
(The  arclength  Is  normalized  so  s“0  at  0«-|-Tr  and 
s  *  1  at  ir/m  -  ^  v.)  The  dashed  curves  are  =  1.39256 
(m  -  3),  1.1^709  (m  -  4),  I.II587  (m  -  5)  and  1.09855  (m  «  6) 
In  Table  III.  The  solid  curves  are  the  limiting  V-states. 

Range  of  existence  of  rotating  V-states  for  2  £  m  £  6.  The 
dots  are  the  lower  end  of  the  range  from  the  incorrect 
analysis  of  Burbea  [ll] . 

Maximum  curvature  versus  R^  for  m  ■  3  rotating  V-states. 
Tl)e  dots  correspond  to  the  third-order  algorithm  (used  with 
■  0.01*  for  R^  <  1.7  and  0.001*  for  R^  >  1.7).  The 
x*s  correspond  to  the  results  In  DO]. 

A  study  of  the  behavior  of  the  m  ■  3  rotating  V-states  near 
R|  using  the  fixed  point  boundary  condition,  (a)  Is  the 

tangent  angle  of  the  contour  at  6  <■  -ir/2.  The  algorithms 
and  discretizations  are:  (A,  dotted  line)  2nd  order  and  1* 
(constant);  (B,  dotted-dashed  line)  2nd  order  and  0.1* 


•  ■ 


Cllnearlj^'  Increasing);  (C,  dashed  line)  2nd  order  and 
0.01®  Cllneariy  increasing) ;  (D,  solid  line)  3rd  order 
and  0.01®  Cline arlj/'*  increasing) .  (b)  The  behavior  of 

selected  contours  near  R*  for  the  3rd  order  algorithm 
and  0.001®  Clineariy' increasing) . 

Figure  11;  A  comparison  of  the  limiting  m  «  3  rotating  V-state  for 

-it/2  £  6  <  -ir/6  using  the  third-order  algorithm  with  0.001® 
(the  dots)  and  the  analytical  formula,  Eq.  (3.23). 

Figure  12;  A  comparison  of  an  analytical  m  ■  3  rotating  V-state, 

“  1.73^93  (the  middle  curve  in  Figure  10b),  with  the 
analytical  formula,  Eq.  (3.30). 


